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Abstract 

We evaluate the decay constants for the B and D mesons and the form factors for 
the semileptonic decays of the B meson to D and D* mesons in a Bethe-Salpeter model. 
From data we extract Vcb = 0.039 ± 0.002 from B D*W and Vcb = 0.037 ± 0.004 
from B Dlv decays. The form factors are then used to obtain non-leptonic decay 
partial widths for B Dir{K) and B DD{Ds) in the factorization approximation. 



1 INTRODUCTION 



In previous papers |]l], ^ we have developed a model for mesons based on the Bethe- 
Salpeter equation (BSE). Recently |0, 3, we calculated the form factors in the semileptonic 
B D{D*)lv decays and extracted the Cabibbo-Kobayashi-Maskawa (CKM) matrix el- 
ement Vcb from data. The key ingredient in the computation of the form factor was the 
construction of the physical states for the B and D mesons in terms of the wavefunction 
obtained by solving a reduced BSE. In this paper we improve upon the work of Ref . in two 
ways. First we establish a theoretical connection between the matrix element of the bare cur- 
rent operator calculated using the model states constructed in and the matrix elements of 
an effective current operator based upon arguments for contributions from neglected config- 
urations. Second, we make an ansatz for the correspondence between the matrix elements of 
the bare and the effective current operator. The effective current operators are then used to 
calculate not only the decay constants and semileptonic form factors similar to our previous 
work but also the branching fractions of non-leptonic decays. 

The discovery of Heavy Quark Symmetry (HQS) in recent years ^ 3 has gen- 

erated considerable interest in the study of systems containing heavy quark(s). It has been 
shown that, in the heavy quark limit, the properties of systems containing a heavy quark 
are greatly simplified. HQS results in relations between non-perturbative quantities, such as 
form factors, for different processes involving transitions of a heavy quark to another quark. 
The development of Heavy Quark Effective Theory (HQET) |^ allows one to systematically 
calculate corrections to the results of the HQS limit in inverse powers of the heavy quark 
mass iriQ. In spite of impressive results obtained in HQET, it has not solved the problem of 
calculating the transition form factors in QCD. In particular, HQS reveals relations between 
form factors but does not provide a determination of the form factors themselves. Further- 
more, the systematic expansion of the form factors in l/mq in HQET involves additional 
non-perturbative matrix elements which are not calculable from first principles. 

We are thus forced to rely on models for the non-perturbative quantities. However, the 
constraints of HQET, which are based on QCD, allow one to construct models which are 
consistent with HQET and hence QCD. We have already demonstrated the consistency of 
our model with the requirements of HQET [Q. 

The parameters in the BSE are fixed by fitting the meson spectrum. Hadronic states 
necessary for the calculation of form factors are constructed with the BS wavefunctions. 
In our formalism the mesons have been considered as composed of qq constituent quarks 
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which defines the hmits of our model space. Our dual thrust in this effort is to correct for 
the limited model space and to carry out new applications. Higher Fock state effects are 
introduced though an ansatz, involving an additional parameter, connecting the bare current 
operator to an effective operator. 

The additional parameter introduced in this process is chosen by a fit to selected exper- 
imental data. An evaluation of the decay constants of the B and D system along with the 
semi-leptonic form factors and the non-leptonic decays is performed without any additional 
free parameter and are, therefore, viewed as predictions of this model. Based on this new 
approach we again extract Vcb from the measured differential decay rate of -B ^ D*W and 
find a 20% increase over our previous results M. We also extract Vcb using recent measure- 



ment of -B ^ Dlu by CLEO [T^. The two resulting values for Vcb presented in this paper 
are consistent with each other. 

The paper is organized as follows: In Section 2 we give a brief review of Bethe-Salpeter 
model for mesons. In Section 3, we discuss the formalism for the calculation of the decay 
constants and the form factors after establishing the connection between the bare current 
oprator and the effective current operator. In Section 4, we discuss non-leptonic decays and 
in Section 5 we present and discuss the results of our work. 



2 BETHE-SALPETER MODEL FOR MESONS 

In Ref. we have developed a model for mesons based on the Bethe-Salpeter equation. 
The wavefunctions for the mesons were solved from three dimensional reductions of BSE, 
called the Quasi-potential equations (QPE). It was found that two reductions give a good 
description of the meson spectrum, including open flavour mesons, over a wide range of 
states. Masses for 47 states were predicted using seven parameters given below with mass 
root mean square deviation of about 50 MeV Ref.p|. 

The interaction kernel in the BSE is written as a sum of a one-gluon exchange interaction 
in the ladder approximation, Vqge, and a phenomenological, long-range linear confinement 
potential, Vcon- In momentum-space this interaction takes the form, 

VoGE + VcoN = -as-, -^ + alini — — -— — - 1 

3 [q — q'Y M^Ocj/i^ —[q — q'y + /x^ 

Here, is the strong coupling, which is weighted by the meson color factor of |, and the 
string tension a is the strength of the confining part of the interaction. We adopt a scalar 
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Lorentz structure Vcon as discussed in 

In our model the strong coupling is assumed to run as in the leading log expression for 

as, 

where j3i = 11 — '^nf/S and Uf is the number of quark flavors, with = M|) ~ 0.12 

where is related to to the meson mass scale through, 

Q' = 7^ML_ + /?^ (3) 

where 7 and (3 are parameters determined by a fit to the meson spectrum. 

In our formulation of BSE there are therefore seven parameters : four masses, m„=mrf, 
rric, mg, mb, the string tension a, and the parameters 7 and (3 used to govern the running of 
the coupling constant. Once the parameters are fixed from the mass spectrum, the meson 
wavefunctions from the BSE can be used to predict physical observables. 

Table 1 shows the values of the parameters used in two reductions of Bethe-Salpeter 
equation referred to as A, B reductions M. 



Table 1: Values of the parameters used in reductions A,B together with root mean square 
deviation from experimental meson masses 





Reduction A 


Reduction B 


nib (GeV) 


4.65 


4.68 


m, (GeV) 


1.37 


1.39 


(GeV) 


0.397 


0.405 


ruu (GeV) 


0.339 


0.346 


a (GeV^) 


0.233 


0.211 


7 


0.616 


0.444 


p (GeV) 


0.198 


0.187 


RMS (MeV) 


43 


50 
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3 Decay Constants and Semi-Leptonic form factors 



The weak decay constants for the heavy hadrons are defined below 

< 0\J^,\P{p) > = ifpPf, 
<0\Jf,\V{p)> = mvfvSf, 



J., 



V„. - A,, 



(4) 



where P and V are pseudo-scalar and vector states and and are the vector and axial 
vector currents. 

The Lagrangian for the semileptonic decays involving the b c transition has the stan- 
dard current-current form after the W boson is integrated out in the effective theory. 

G p 



w 



2^2 



l46C7/.(l-75)&^>7'^(l-75)/ 



(5) 



The leptonic current in the effective interaction is completely known and the matrix element 
of the vector (V^) and the axial vector {A^) hadronic currents between the meson states are 



represented in terms of form factors which are defined in the equations below |[TT| 



<D{pd)\J^B{pb) > 



{Pb+Pd)^ - 



+ 



2 2 

m|j — m 



D 



(6) 



where q = Pb — Pd- 



<D*{p')\J^\B{p)> = h,e^^''^elp^p'p + he*^ + h2{p + p'r + Hky 



(7) 



with 



hi 



2V{k^) 
rriB + tud* 
i{mB + m£).)y4i(A;^) 

TUB + niD* 
2mD4Ao{k') - A,{k')) 



ie* ■ k- 



A;2 



{mB + 'm£)*)Ai{k'^) — {niB — m£,*)A2{k'^ 



2m 



D* 



5 



where k = pb — Pd*- Fq, Fi, V, Aq, Ai, A2, and ^3 are Lorentz invariant form factors which 
are scalar functions of the momentum transfer (Pg — P£){Pd*)Y ■ The calculation of the 
decay constants and form factors proceeds in two steps. In the first step, the full current 
from QCD is matched to the current in the effective theory (HQET) at the heavy quark mass 
scale [0 . Renormalization group equations are then used to run down to a low energy scale 
yU ~ 1 GeV where the constraints of HQET operate and where it is reasonable to calculate 
matrix elements in a valence constituent quark model like the one we employ here |jT^. We 
have already described the first step in our previous publication and therefore we will not 
repeat it here. 

The second step is the calculation of the matrix elements of the currents in the model to 
obtain the decay constants and form factors . Such a calculation requires the knowledge of 
the meson wavefunctions. In our formalism the mesons are taken as bound states of a quark 
and an antiquark. The wavefunctions for the mesons, as already mentioned, are calculated 
by solving the Bethe-Salpeter equation §]. We construct the meson states as 



|M(Pm, ^,"^J)) = ^/2M^ J £p{LmLSms\Jmj) {smsSrUslSms) 



Tfl- Tfl 

*Lm,(p)|g(^PM - P, m,-))|g(^PM + P, m,)) (8) 



where 



M 
En 



\ 



{Eg + niq 
2m„ 



X 



{Eq+niq) 



rriq + TUg 



Iml + p2 (9) 
and Mh is the meson mass. The meson and the constituent quark states are normalized as 
(M(P'M,J',mS)|M(PM,J,mj)) = 2E5\^'^ ^jd^'mj (10) 



(g(p',m'J|g(p,m,)) 



^(5='(p'-p)5„,, 



'111 



ma 



In constructing the meson states we maintain a constituent quark model approach as 
we do not include qq sea quark states nor the explicit gluonic degrees of freedom. We also 
assume the validity of the weak binding approximation |TB|, M. In the weak binding limit 



our meson state forms a representation of the Lorentz group, as discussed in Ref . |I3| , if the 
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quark momenta are small compared to their masses. Assuming that the quark fields in the 
current create and annihilate the constituent quark states appearing in the meson state, the 
calculation of the matrix element of the current operator then reduces to the calculation of 
a free quark matrix element. In the rest frame of the B meson with a suitable choice of the 
four- vector indices in Eqs.(6,7) we can construct six independent equations which we can 
solve to extract the six form factors. 

This model space representation may be viewed as the leading characterization in an 
expanded representation which more accurately represents the exact states. We assume that 
the effects of Higher Fock states, representing gluons or sea quarks, in the calculation of the 
matrix element of the bare current operator are represented by the matrix element of an 
effective operator in the model space. In other words, with the notation "e" labelling exact 
states, and "m" labelling model states. 



\rJf\Mr{Pi) > 



'12) 



where the higher Fock state effects are included by the following replacement in the calcula- 
tion of the matrix element 



(13) 



In the above, p, p' are the internal momenta of the quarks in the initial and final mesons 
and $i(p), $2(p') are the initial and final meson wavefunctions. We will use the very simple 
ansatz 



n{p) 



(14) 



We will fix the parameter a by simply fitting to the available experimental data and 
lattice results of the decay constants of the leptonic decays. Note that we will use the 
same value of a for decays involving B and D decays. This is consistent with Heavy Quark 
Symmetry. 

The expressions of the decay constants in terms of the wavefunctions are given as |TB 



fi 

Fp{p) 

Mp) 



12 r°° p^dp 
M Jo '2^\ 

p2 



iTTlr, 



p" 



[niq + Eq){mq + Eg] 

AEqEq 



i)p{p) 

tjjvip) 



F^{p) 



Eg){mg+Eg 



3{mg + Eq){mq + Eq 



(15) 
(16) 
(17) 
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where ipp{v) are the wavef unctions of the exact states. Using Eq.(12), we can then obtain 
the form factor in terms of the BSE wavef unctions. 



4 Non-leptonic Decays 

Non-leptonic decays arise from W exchange diagrams at tree level. Strong interactions 
play an important role in these decays by modifying the weak vertices through hard gluon 
corrections and then the long distance QCD interactions result in the binding of the quarks in 
the hadrons. An effective Hamiltonian of four quark operators is constructed by integrating 
the W-boson and the top quark from the theory. The effects of the short distance and the 
long distance QCD interactions are separated using the operator product expansion where 
the Wilson coefficients account for the short distance effects while the long distance effects are 
incorporated in the matrix element of the four quark operators. The effective Hamiltonian 
operator for 6 — c transition can be written as 

Heff = ^V.^V:,[cM0, + cM02] 

01 = da^{l - 75)MiCj7''(l - 75)6j 

02 = da^.{l-l5)ujCj-f^{l--f5)bi (18) 

where i and j are the color indices. The Wilson's coefficients Ci and C2 at the scale /j, = mi, 
have values 1.132 and -0.286 respectively [^. 



The matrix element of a two body leptonic decay of the type B XY requires the 
evaluation of the matrix element 

M = <X,Y\Heff\B> 

where -ffe// has a current x current structure. The matrix element is usually calculated 
using the factorization assumption where one separates out the current in Hgff by inserting 
the vacuum state and neglecting any QCD interactions between the currents. The matrix 
element above written as a product of two current matrix elements is 

M ~ < X|J^|0 >< r|J'^|5 > (19) 

In B decays, for e.g B D^n^ , the energetic quark- ant iquark pair in the pion is created 
at short distance and by the time it hadronizes it is far from the other quarks so it should 
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be a good approximation to neglect the QCD interaction between the two currents creating 
the final sate particles. A detailed description about the validity and the corrections to the 
factorization assumption can be found in Ref.|p.(j||. 

In this paper we will look at decays where the particle Y is a. D oi a. D* meson because 
one can then use the semi-leptonic form factors calculated in the previous section to compute 
< y| J'^|i? >. The X will be either a light meson( tt,K,p,K*) or a D{D*) meson. For the 
light mesons the decay constants fx =< X\J^\() > are available from experiment while, for 
the heavy mesons, we will use the decay constants calculated in the previous section. 

The expressions for the square of the matrix element for the processes D[D*)'k[p) 



are 



,2^2, 



V2 



— m 



2 \2 

d) 



where Nc is the number of colors. 

|M|2(5°^DV) - 



where p is the momentum of the decay products in the rest frame of the B. 



and finally, 

|M|2(5° ^ D*+p- 



,G 



,2\2 2 P 



(20) 



(21) 



(22) 



D* 



T3 
T4 



8V^ 



{ttlb + rriD' 



2 2 

-p TUb 



2 + 



{irLBEo* 



m 



D*) 



4 4 

p 



AA1A2 

Similar expressions can also be written down for the B DD decays. 



El 



El, 



mi m 



ED*Ep{mBED* 



m 



(23) 



As in Ref. |^B[ we will include the effect of corrections to the factorization assumption by 
the replacement 



ai 



ai 



Ci(/i) + 



C2(/^) 



+ C2{p)e^{p) 
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Table 2: Decay constants of the B and D mesons in MeV 



Decay Constants 


Our Results 


Lattice Results [T^ 


fo 


209 


196(9)(14)(8) 




237 






213 


211(7)(25)(11) 




242 




fB 


155 


166(11)(28)(14) 




164 





The nonfactorizable corrections and s^{^) are defined in Ref.|T^ and may be process 

dependent. We will, however, treat ai as a process independent free parameter that we will 
fit to data. RGE analysis suggests that oi ~ 1 + Q{l/Nl). 



5 Results and Discussions 

In previous papers 0, ^, a covariant reduction of the Bethe -Salpeter equation (BSE) was 
used to calculate the Isgur-Wise function. The BSE was solved numerically and the param- 
eters appearing in it (the quark masses, string tension and the running coupling strength for 
the one gluon exchange) were determined by fitting the calculated spectrum to the observed 
masses of more than 40 mesons. The resulting mass spectrum of the analysis was found to 
agree very well with the experimental data. Once the parameters of the model were fixed, 
the meson wavefunction could be calculated from the BSE. This wavefunction was used to 
calculate the Isgur-Wise function and determine Vcb [§]. 

In our present approach we evaluate the decay constants, the form factors for the semilep- 
tonic decays B —>■ D*W and B —>■ DW with the effective current operator defined in Eq.(13) 
treating a of Eq.(14) as an adjustable parameter. The value of a is fixed by fitting the 
leptonic decay constants. We find a = 0.7GeV^^ provides a good fit and use this value in 
all the calculations in this paper. 

We present our results for the decay constants of the heavy mesons in Table 2. For the 
sake of comparison we also show lattice calculations of the decay constants. The errors in 
the second column of the Table are, respectively, (1) the statistical errors; (2) the systematic 
errors of changing fitting ranges, as well as other errors within the quenched approximations; 
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and (3) the quenching error. The results in Table 2 show that our calculated decay constants 
are similar to the lattice results. 

On the other hand, our calculation for the decay constants of the light mesons tt, K etc 
are not in good agreement with the experimental numbers. In fact, the light meson decay 
constants are larger than experiment by a factor of 2. This is not surprising as our formalism 
is designed for the heavy meson system. 

In Fig.l we show the form factors Fq, Fi, V, Aq, v4i, and A2 as a function of q^. In Fig. 2 
we show a plot of the differential decay rate for B — > D*W. We obtain a good agreement 
with the shape of the experimental data and extract \Vcb\ = 0.039±0.002. This is within 
the range of the presently accepted values for \Vcb\ [p^ ■ 



For the decay B — > DW, in Fig.3 we show a plot of F(a;)|K;fe| versus u where the data 
points are taken from measurements reported in Ref.[|lO]. The variable u = (M^ + M|, — 
q^) / {2MbM£,) where is the invariant mass squared of the lepton neutrino system. We 
extract \Vcb\ = 0.037 ± 0.004 by a fit to the data in Fig.3. 

Note that the values of Vch extracted from the two different experiments are consistent 
with each other. As a further test of our formalism we present our calculations of the 
nonleptonic decays of the B meson to DD and DK{7c) final states. Experimental values 
of some of the decays are already available and new results are expected soon. We present 
our results in Table 3. The values of the light decay constants used in our calculations are 

= 130MeV, fx = 159 MeV, fx* = 214 MeV and fp = 208 MeV. 

The parameter ai calculated on the basis of a fit has the value 0.88 which is close to 
1 as is expected from RGE analysis which gives ai ~ 1 + 0{1/N^) suggesting a value for ai 
in the range 0.9 — 1.1 

/^From Table 3 we find a good agreement of our calculation with data for the D*Ti{p) 
final states. Our results for the DK final states are quite similar to those in Ref. |p.6||. 
This continues to be true for the DD and D^D~ final states. Our results for D^D~ and 
D^D~ final states are somewhat smaller than the central values from experiment though the 
measurements have large errors. 

Combining the experimental errors in quadrature the difference between theory and ex- 
periment is less than 1.5a in all cases but the theory predictions are systematically lower 
for these cases. It appears that as we increase the mass of the decay products, as in the 
DD final states, and decrease their kinetic energy the expected deterioration of the factor- 
ization approximation may be showing up through a systematic difference between theory 
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Table 3: Non-Leptonic Decay Rates for B meson 



r 1 UL^COO 


V 7 ILL rvcoUlLo 


LCL^li IN C UUci L 1 -L VJ 1 


J— /A.JJ L |ZiU|| 




yJ.O'-tO 


\J .OKJyJ 


31 Oro 040V0 090"! 




yJ.OO ± 


990 


980l'0 040V0 01 a) 


J_J 7 ±y U 


n 7QQ 


750 


840l'0 1 60 VO 070"! 


U 7 ±y u 




850 


730l'0 1 50U0 030"! 


RO . n+ r^-- 

1 J 7 J _y J V 




n 90 




J J 7 J _y j\ 




90 




RO . n+K*^ 

J J 7 J _y J V 




40 




RO . n*+ 

J J 7 J _y ± \ 


n 4Q 


50 








40 




RO . n*+ n- 


n 99 


"^0 




RO . n+ n*- 

1 J 7 I J I J 


n 97 


"^0 




DO . n*+ n*- 

i J 7 I J I J 




RO 








1 O'^O 


740l'0 99VO 1 R\ 


RO . n*+n- 


n 490 


700 


Q4l'0 94VO 9'^'l 


RO ^ n+n*- 




950 


1 140f0 42V0 2H) 


50 ^ D*+D*- 


1.35 


2.450 


2.0(0.54)(0.05) 


B- D^D- 


0.33 


0.40 




B- D*^D- 


0.210 


0.30 




B- D°D*- 


0.27 


0.40 




B- D*^D*- 


0.64 


0.90 




B- -> D^D; 


0.829 


1.090 


1.360(0.280)(0.330) 


B- D*^D- 


0.552 


0.750 


0.940(0.310(0.23) 


B- D^D*- 


0.696 


1.020 


1.180(0.36)(0.29) 


B- D*°D*- 


1.830 


2.610 


2.700(0.810)(0.660) 
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and experiment. This motivates future efforts to examine corrections to the factorization ap- 



proximation [0. We have resisted the temptation to allow ai to have a process dependence 
even though two values for ai would yield an excellent description of the known nonleptonic 
decay rates. It is trivial to relax this restriction if the reader chooses to do so. 

In conclusion, we have presented the calculation of form factors and differential decay 
rates in B D{D*)lu transitions in a Bethe-Salpeter model for mesons. The parameters 
of the bound state model were fixed from the spectroscopy of the hadrons. The effects of 
higher Fock states in the hadron state were included in the definition of effective current 
operators. A simple ansatz connecting the effective current operator to the actual current 
operator was used involving only one parameter. After adjusting this parameter to fit certain 
decay constants, we found good agreement with data and extracted \Vcb\ = 0.039 ± 0.002 
from B D*W and Vet = 0.037 ± 0.004 from B Dlv decays. Calculations of the decay 
constants of the B and D mesons were also performed with results that are similar to lattice 
results. Finally, the form factors were used to evaluate the non-leptonic B D7[{K) and 
B — s> DD{Ds) decays in the factorization approximation and good agreement was obtained 
with data. 
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6 Figure Captions 

Fig.l: The calculated form factors Fo,Fi, V, Aq,Ai, and A2 as a function of g^. 

Fig. 2: The differential decay rate for B D*W with and without the QCD correction, 
together with the corresponding values of Vet- Data from Ref.|0. 



Fig. 3: F(ti;)|V^b| versus uj for B ^ Dlu. Data from Ref . [p!Of 
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